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Abstract 

The effect of point defects on persistent currents in mesoscopic 
rings is studied in a simple tight-binding model. Using an analogy 
with the treatment of the critical quantum Ising chain with defects, 
conformal invariance techniques are employed to relate the persistent 
current amplitude to the Hamiltonian spectrum just above the Fermi 
energy. Prom this, the dependence of the current on the magnetic flux 
is found exactly for a ring with one or two point defects. 

The effect of an aperiodic modulation of the ring, generated through 
a binary substitution sequence, on the persistent current is also stud- 
ied. The flux-dependence of the current is found to vary remarkably 
between the Fibonacci and the Thue-Morse sequences. 
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The study of the critical behaviour of spin systems close to surfaces and 
interfaces has revealed a very rich phenomenology and has provided much in- 
sight into new phenomena. Two-dimensional systems have been of particular 
theoretical interest, because of fruitful connections with conformal invariance 
techniques, see fl], ||, |^, Q for an introduction and review. Here, we present a 
new application of some conformal invariance methods to a simple problem 
of persistent currents, which we now specify. 

Consider a small ring, made from a semiconducting or a normally con- 
ducting material. If the ring is threaded by a magnetic flux, one observes 
experimentally a persistent current, provided the ring has mesoscopic dimen- 
sions, that is, its diameter is in the nanometer range. Measured examples 
include Au or GaAS-AlGaAs rings. For larger rings, the current de- 
creases and is no longer perceptible for rings of macroscopic size. The problem 
has received much attention. Non-interacting rings were studied in the classic 
papers by Kulik [^] and Biittiker et al. 0] . Disordered rings were considered 
in ip. On the other hand, persistent currents for interacting ID electrons 
were first investigated in |0|, 0]. See also p|, |l3Hl|, |15|, ||, [l^, |18], [H and 



references therein. We refer to for a review and to the nice book by Imry 
for a detailed presentation of the physical background. 
Here, we shall limit ourselves to a simple toy model which describes non- 
interacting (spin-less) electrons moving around a small, thin ring penetrated 
by an Aharonov-Bohm flux $, see Figure ^ The flux $ = / Adx = AL, 
where A is the vector potential and L is the length of the ring. Assume that 
the ring is discretized into cells of size a. We use a simple tight-binding 
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Figure 1: Aharonov-Bohm geometry. The places of the defects and their 
couphngs Ki 2 are indicated. The magnetic field is perpendicular to the plane. 



model which describes the hopping of electrons from one cell to the next 

n=l L \ 9oLJ J 

where thop is the hopping rate and $o = hc/e is the flux quantum. Hamilto- 
nians of this kind have been proposed earlier, see [jTT], [1^, to study the effect 
of an Aharonov-Bohm flux on a Luttinger liquid. Our model (|l|) is, with 



respect to [0, |T^, simplified even further in neglecting the charge- density 
interactions and only considering the zero-temperature case. On the other 
hand, we allow for point-like defects to occur in the ring, which phenomeno- 
logically lead to a local modification of the hopping rate, as already studied 
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We shall limit ourselves throughout to the case of half-filling, when exactly 
M = N/2 particles are present and assume that is even. The flux $ can 
be related to the phase 

= 27r$/$o (2) 



of the wave function through the gauge transformation 

^jJ{n) — exp (^^^j^''^^ 'ipin) (3) 
leading to the conventional tight-binding Hamiltonian 

N 

H = -thop E [^^(^ + l)^(^) + ^^(^)^(^ + 1)] (4) 

n=l 

with the boundary condition ip{N+l) = e'^'^ipi'^) ^i-nd where is the number 
of sites. This is nothing but the classical Byers-Yang-Bloch result that the 



system is periodic in the flux $ with period $0^ see |TT|, [19], |20|, ^ and 
references therein. In particular, the amplitude j of the persistent current is 
obtained from the ground state energy Eq{^) of H as 

It is well known that for a perfect ring, j = (thop/7r)(0/A^) inside the period 
interval — vr < </> < tt and continued periodically outside it ||20|, Note 



the scaling of j with the ring length L = Na 0) so that the current 
vanishes in the bulk limit L cxo, as expected. 

We now consider the situation where in an otherwise perfect ring there 
is a single link where the hopping rate is thopP^Q rather than thop, see Fig- 
ure |l|. By rotational symmetry, we thus have a defect boundary condition 
iIj{N -|- 1) = with a complex defect strength k, = pe"*"^. Formally, the 

diagonalization of the free fermion XX Hamiltonian (0) proceeds along ex- 
actly the same lines as had been developed earlier for the treatment of the 

■^To facilitate comparison between the models considered here, signs and phases are 
arranged such as to reproduce this result for the defectless case throughout. 
"^The impurities studied here are only forward scattering. 
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quantum Ising chain in the presence of point defects (or, equivalently, the 
classical 2D Ising model with line defects) Through a Jordan- 



Wigner transformation and a subsequent canonical transformation, H can 
be written in terms of free fermionic harmonic oscillators rjk 

^ = EA.(r]k-^) (6) 

where = — thop cos k and the values of k are determined by solving 

p2 sin(A;(A^ - 1)) + 2p cos(/)sin A; — sin(fc(A^ + 1)) = (7) 

The physical ground state is obtained by filling up the lowest M = N/2 
one-fermion states. Formally, the ground state energy — can now be 

calculated and the current j be obtained. Practically, this is not straightfor- 
ward, since there is no closed solution to for |p| ^ 1. 

However, to a collective effect as a persistent current only the lowest ex- 
cited levels should make a significant contribution. Plotting A^ as a function 
of k, one sees that for half-filling, the lowest modes occur at values of k 
around | and Linearizing A^ close to the Fermi energy, we can therefore 
replace the exact H by an effective low-energy Hamiltonian ifgff of the form 



6N 

(8) 

where n^^^ are fermionic number operators representing the two lowest modes 
and where the one-fermion energy shift isV = T>{p, 0), such that T>{1, 0) = 0. 
The same form had been found before for the Ising quantum chain with 
defects ||, |3 . 



4 



In (^) we have dropped a constant contribution of the form FN + S. 
However, these terms have no bearing on the persistent current. To see 
this we use the analogy of the XX chain (^) with the Ising quantum chain. 
In turn, the Ising quantum chain can be viewed as the Hamiltonian limit 
of a two-dimensional classical Ising system, making contact with standard 
thermodynamics. From this, we recognize F as the bulk free energy density 
per site and S corresponds to a surface free energy term. This interpretation 
is carried over to the XX chain. Since the bulk free energy density F is 
independent of any boundary condition, it is thus also independent on the 
complex defect strength k. On the other hand, in the same framework the 
surface free energy S may be sensitive to boundary effects. Because of the 
gauge transformation the flux can also be viewed as a bulk interaction, 
and it follows that S will merely depend on the real defect strength p, but not 
on 0. In conclusion, since F, S are both independent of 0, eq. contains 
all the information required for the determination of the persistent current 
through (H). 

Because the dispersion relation close to the Fermi energies is linear, the 
effective Hamiltonian represents a conformally invariant system p, |^, ^ . Con- 
formal invariance allows to relate the (half- filled) physical ground state energy 

Eo= Afc = 2^o,cfr (9) 

to the constant c^s and furthermore to the one-fermion energy shift "D. Here, 
Eo^eS is the ground state energy of H^s and the factor 2 comes from the 
symmetry of at half-filling. 

One of the main results of 2D conformal invariance (e.g. 0, 0, Q), applied 
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to critical classical spin sytems, states that it acts as a dynamical symmetry 
for the associated quantum chain of N sites, with Hamiltonian 

where Lq, Lq are generators of the Virasoro algebra and c is the central 
charge. Here, the quantum Hamiltonian H is obtained from the transfer 
matrix T = exp(— rif) of the classical 2D spin model, see ^, As it 
stands, (p^Of ) applies to systems without any defects at all and in particular to 
quantum chains with periodic boundary conditions. One may easily recover 
the spectrum-generating conformal algebra for the excitations described by 
i^eff above the Fermi threshold, in complete analogy with the familiar con- 
struction for critical Ising spin systems, see 

Now, introducing non-periodic boundary conditions (which may also be 
viewed as a point defect in an otherwise periodic quantum chain), it can 
be shown that the transformations generated by the diagonal subalgebra 
Kn = Ln + Ln are preserved while the other generators L„ — L„ no longer 
correspond to a symmetry of the system. Such a point defect in the quantum 
chain Hamiltonian corresponds, via a conformal transformation, to a semi- 
infinite line defect in the associated 2D spin model, see ^ and references 
therein. Specifically, for the Ising quantum chain it turns out that the per- 
fect generators L„ are modified into Virasoro generators Ln{T>) depending 
explicitly on the one-fermion shifts such that [E^ 

= K^V) = L„(P) + + ^VX,o (11) 
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satisfy a Virasoro algebra with central charge c = 1 and 

= I {k,(V) - ij) (12) 

which generalizes (|T0|). Although invariance under the full conformal algebra 
generated by the set L„} is broken, the diagonal conformal subalgebra 
generated by the Kn still acts as a dynamical symmetry. Thus, Ces = \ ~^'^^ 
[ p^ and the only change in the ground state energy of H which can give rise 
to a persistent current is 

5Eo = 25Eo,eff = (27r/iV)thop^^2 (13) 

This is the desired result, which relates the current j to the shift T> of the 
lowest excitations. 

To calculate V to leading order in we first shift k to A; + 7r/2 in 

order to have = at the Fermi surface and then lei k = a/N in 
This gives cosa = (— l)^'^[2p/(l + p^)]cos0. From (|]), we also have k ~ 
{2TT/N){i +1/2±V), with £ = 0, 1, 2, . . .. Using these two formulas for k, 
we find V = T>{p, 0). From (^) the current j = —d{5EQ) / d(j) is obtained as a 
function of the phase 0. Explicitly, 

., u arccosfu cos</)) sin(/) 1 2p . 

= thop ^=4=J^- , u = 14 

TTV 1 — M COS^ (p 1\ I + p 

and we find again the scaling j ~ A^~^. Note that this expression is invariant 
under the transformation p ^ 1/p which reflects that the defect, even for 
p > 1, disturbs the coherence of the many-particle system. As it stands, (0) 
is correct for M odd. For M even, merely replace (p ^ (f) — n. The size- 
independent current amplitude j{4>)N is shown in Figure |^a for an entire 
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Figure 2: Persistent current amplitude j{(p)N in the tight-binding model for 
(a) a single defect with defect strength k, — pe~"^ for M — N/2 odd. For 
two opposite defects with strengths ki = pe~^'^ and H2 = P the persistent 
current is shown for (b) N/2 even and (c) N/2 odd. The curves correspond 
to different values of p as indicated and thop = 1- 
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period of the phase and for several values of the defect coupling p. While 
the periodicity with is kept unchanged with respect to the perfect case, 
the amplitude decreases rapidly with increasing values of |p — 1|. Of course, 
this result can be checked through direct numerical diagonalization of the 
tight-binding Hamiltonian (^) and may also be reproduced from a simple 
continuum description of the model . 

For two defects which are placed opposite to each other (Figure |I]), we 
find in an analogous fashion 

a arccos (a cos ± &) sin 1 



j±(0) = 



op 



TT^l - (a COS 0± 6)2 N 

4piP2 , {1-pI)(1-pI) 



(15) 



(i + P?)(i + pi) ' (i + pf)(i + Pi) 

where + corresponds to N/2 even, — corresponds to N/2 odd and pi^2 are 
the couplings at the defects. For pi = 1 or p2 = 1, one simply recovers (|1^). 
In Figure ^dc, we show the current amplitude for the number of particles 
M = N/2 even and odd, respectively and for pi = p2 = p- In contrast 
to the single-defect case, there is a clear parity effect. For an even number 
of particles, the amplitude is broadly similar to the single-defect case (up 
to a renormalization of p). In particular, the current j(0) is a continuous 
function of over the entire period. For an odd number of particles, the 
flux-dependence of the current is almost linear for a wide range of values of p 
and there is a jump in j at the ends of the period interval. Also, the presence 
of the defects merely changes the prefactor. 

We now consider the case of two defects placed anywhere on the ring. 
The first one may always considered to be placed at the origin, while the 
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second one will sit at the site R = XN. The continuum dispersion relation 
remains unchanged, but (0) has to be replaced by 

(l - pj) (l - pI) cos [a{l - 2A)]±(l + pj) (l + pI) cosa±4pip2(-l)^^ cos0 : 

(16) 

where = 2M is the number of sites, k = a/N (after the shift k k + TT/2) 
and +(— ) corresponds to R even (odd). Since A = R/N is rational, we 
may write X = p/q with p,q mutually prime. Eq. ([TBD then becomes a 
polynomial equation of degree 7(g) in the variable cos(2a/g) = cos(2'7rD'-*)), 
where 7(g) = g/2 if g is even and 7(g) = g if g is odd. The current is then 
written in terms of the 7(g) distinct one-fermion shifts as 
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It is straightforward to extend the treatment to an arbitrary number of point 
defects, following the lines of |2^ for the Ising quantum chain. 

These results were derived by going back to the explicit diagonalization of 
the underlying tight-binding model and using the analogy with the quantum 
Ising chain for defect boundary conditions. In the Ising quantum chain, a 
defect boundary condition corresponds to a marginal perturbation |^ which 
leads to a shift V in the one-fermion energies which depends continuously 
on the defect strengths and consequently to a continuously varying effective 
central charge Ces- Extending existing field-theoretical treatments of defect 
lines in 2D models |2^, |2B|, to complex defect strengths might provide 

further insight into the problem at hand.0 

^For example, the generators (p3| ) can also be obtained using the automorphisms of the 
underlying Kac-Moody algebra p8| . This more abstract construction is independent of a 
free-fermion realization of the model. 
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Having seen that at least in the tight-binding model with point defects, 
the persistent currents are related to the shift in the low-lying excitations 
above the Fermi energy, we now briefly discuss persistent currents in ape- 
riodically modulated rings threaded by a magnetic flux. The Hamiltonian 
is 

N 

H = -Y.tn [i^\n + l)i;{n) + ^p\n)^{n + 1)] (18) 

n=l 

at half-filling, together with the flux boundary condition used before and 
where the aperiodic modulation in the hopping matrix element t„ is generated 
through a binary substitution sequence p9|, |30|. We shall consider here 



Fibonacci : A ^ AB , B ^ A 
Thue-Morse : A ^ AB , B ^ BA (19) 

and the sequences are obtained by starting from the letter A and iterat- 
ing the above procedure many times. For example, using the Thue-Morse 
substitution rules, the first three iterations give the following sequences 

A 



i = 1 
i = 2 
i = 3 



AB 
ABBA 
ABBABAAB 



(20) 



A modulation in (0) is obtained by letting t„ = thop if the site n is the letter 
A and t„ = thopP otherwise. For such aperiodic modulation the cumulated 
deviation from the average hopping rate t, A(L) := I]fc=i(^fc — ^) scales 
with the chain length L as A(L) ~ , where uj is the so-called wandering 
exponent, caracteristic of the sequence |^ If < 0, as it is the case for 



the two sequences under consideration in this paper, A(L) remains bounded 
for L — > oo. 
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Figure 3: Current amplitude Nj for a ring modulated with the Fibonacci 
sequence, N = 144, thop = 1 and different values of p. 



For critical spin systems. Luck |3T| has formulated a criterion which for a 
given spin system and a given sequence predicts whether the aperiodic modu- 
lation is relevant, marginal or irrelevant. According to the Luck criterion, the 
modulations generated by both the Fibonacci and the Thue-Morse sequence 
are irrelevant for the Ising quantum chain and the low-lying excitations of 
H are in the N oo limit not affected by the modulation. Because of the 
close connection of the tight-binding model with the critical point of the Ising 
quantum chain, we expect to retain the scaling j ~ A^~^ of the persistent 
current, but we are curious about the flux- dependence of 

In Figure 0, we show the current amplitude jN as obtained for several 
values of p for an aperiodic ring modulated with the Fibonacci sequence. 
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Figure 4: Current amplitude Nj for a ring modulated with the Fibonacci 
sequence, p = 1.5 and varying number of sites A^. 

We also include the straight line behaviour for the unmodulated ring for 
comparison. We observe that the modulation reduces the current density 
with respect to the unmodulated case and that the dependence of j(0) on 
the flux is similar to the one found for the single defect in Figure |^a. In 
Figure ^, we illustrate for fixed p the convergence of the finite-size data with 
A^. We see that quite large lattices are needed to achieve convergence of 
the current amplitude. Finally, in Figure |^ the current amplitudes for the 
Thue-Morse sequence are displayed. In this case, the convergence with N is 
much faster than for the Fibonacci sequence and already for = 64, there 
is no perceptible change of the amplitude with N. 

Remarkably, the current depends on the flux in very different ways for the 
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Figure 5: Current amplitude Nj for a ring modulated with the Thue-Morse 
sequence, N = 64, t^op = 1 and several values of p. 

two sequences, although we observed indeed the expected scaling j ~ A^~^. 
For the Fibonacci sequence, the amplitude Nj{4>) with p ^ lis always smaller 
than the amplitude for the non- modulated ring p — 1, independently of 
whether p is larger or smaller than unity. That is not so for the Thue-Morse 
sequence. Rather, the amplitude Nj is weakened for p < 1, but for p > 1, 
the amplitude becomes considerably larger than in the unmodulated case. In 
fact, for p oo, we recognize a limit form Nj{(j)) = — sin (p. In spite of being 
irrelevant to the critical behaviour of Ising quantum chains, the details of the 
aperiodic modulation are important for the flux-dependence of the persistent 
current. These finer properties are not reflected in the low- lying excitation 
spectrum of the Hamiltonian. 
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In conclusion, the role of defects on the persistent current in mesoscopic 
rings was studied through a simple tight-binding model. From the analogy 
with the local critical behaviour of conformally invariant spin systems with 
defect lines, when extended to complex defect strengths, explicit formulae 
for the persistent current were obtained. The fiux-dependence of the current 
shows an unexpected parity effect. In addition, the persistent current was 
found for two types of aperiodically modulated rings. The flux-dependence 
of the current cannot be obtained from a consideration of the low-lying 
excitation spectrum alone, in contrast to the case of point defects. It would 
be interesting to sec to what extent these observations can be generalized 
beyond the simple tight-binding model. 

Acknowledgement: It is a pleasure to thank H. Schoeller for very useful 
discussions. 
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